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Abstract 

We use the holomorphic anomaly equation to solve the gravitational correc- 
tions to Seiberg-Witten theory and a two-cut matrix model, which is related 
by the Dijkgraaf-Vafa conjecture to the topological l?-model on a local Calabi- 
Yau manifold. In both cases we construct propagators that give a recursive 
solution in the genus modulo a holomorphic ambiguity. In the case of Seiberg- 
Witten theory the gravitational corrections can be expressed in closed form as 
quasimodular functions of T(2). In the matrix model we fix the holomorphic 
ambiguity up to genus two. The latter result establishes the Dijkgraaf-Vafa 
conjecture at that genus and yields a new method for solving the matrix model 
at fixed genus in closed form in terms of generalized hypergeometric functions. 
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1 Introduction 

The holomorphic anomaly equations, discovered inpQ in a world-sheet analysis of 
a topological twisted u-model known as B-model, are a generalisation of Quillens 
anomaly to higher genus and more general world-sheet theories. Gauge theory are 
embedded in string theory and in the N = 2 context the latter can be obtained 
in double scaling decoupling limit. The holomorphic anomaly equations commute 
with that decoupling limit and the recursive procedure which determine the higher 
PP, which describe certain terms in the coupling to gravity, can be build 
up entirely from gauge theory quantities. In section the corresponding topological 
partition of the 4d N = 2 SUSY gauge is determined by solving the holomorphic 
anomaly recursively, up to finitely many terms, which have to be fixed by analyzing 
the boundary behaviour of the F^ g \r,f). The properties of F^ 9 \t,t) are to a large 
extend fixed by the modular group of the corresponding Seiberg-Witten curve and 
we are able to write global expressions them in terms of "almost holomorphic" mod- 
ular functions of this group. Various holomorphic limits are readily taken from our 
expressions and provide conjectural solutions to the unitary matrix model [3]. 
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Riemann surfaces £ 9 of all genus can be embedded in noncompact Calabi-Yau 
threefolds X, i.e. the local limit of the 5-model exist for an arbitrary Riemann surface. 
There are choices in this embedding, which affect the reduction of the holomorphic 
(3, 0) form from X to one form A on S 9 , which is a key datum of the local .B-model 
limit. Except for the product case, A will be a meromorphic form with residua and 
eventually boundary data on open T, 9t h- In the pure SW-case, there are poles with 
vanishing residua. In the massive SW-case and other local geometries, one one would 
have parameters associated with the non- vanishing residua. In section |3] we apply the 
holomorphic anomaly equations to iV = 1 4d gauge theory in a geometry proposed 
by j3]. This is a case where A has essential singularities or differently put one deals 
with open Riemann surfaces. According to the Dijkgraaf-Vafa correspondence [H| the 
holomorphic anomaly equations provide here in a particular region solutions to the 
large N expansion of a complex matrix model, which is the solution to an open string 
problem. This has a description in the A and the B language as summarized following 
commuting diagram of duality relations 



open topological string with 

3d Lagrangian cycles 

-> 3d Chern Simons theory 



large N-duality 



geometric transition 



closed topological A-model 
on blown up geometry with 
even form R flux 



mirror 
duality 



mirror 
duality 



holomorphic Chern-Simons 
th. reduced on holom. curves 
> matrix model 



large N- duality 



geometric transition 



closed topological B-model 
on deformed geometry with 
odd form R flux 



The Dijkgraaf-Vafa conjecture, relevant to calculate effective terms in AT = 1 super- 
symmetric theories in 4d, involves non-compact Calabi-Yau geometries, which are in 
particular not toric. Explicit test of the conjecture are therefore difficult and have 
been done only at tree level [Sj and at genus one [HI CHj • In the genus one case it is 
found that the free energy for multi-cut matrix models can be written in closed form 
[HUH]. We set the B-model formalism up to get recursively closed expressions at all 
genus and check the Dijkgraaf-Vafa correspondence of the first non-trivial example, 
the two cut case, explicitly at genus two. 

As in the N = 2 case, we make a detailed analysis of the moduli space and the 
modular transformation of the periods in Appendix [21 which enables us to solve the 
theory at various regions in the moduli space. Different from in the N = 2 case 
however there is a divisor with essential singularities of the periods in the moduli 
space, where the perturbative description breaks down. From the topological string 
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point of view we obtain solutions for a new class of non-toric, non-compact Calabi-Yau 
geometries, which have no point of maximal unipotent monodromy. 



2 Gravitational corrections to N=2 Seiberg-Witten 
Gauge theory 

In this section we introduce a B-model, which calculates the higher genus space-time 
instantons of N = 2 Seiberg-Witten gauge theory. The genus g generating function 
Fb)( a ) f these space-time instantons describe gravitational couplings to the gauge 
theory ^3]. The coefficients of 

oo 

F = J2x 29 - 2 J ri9) (a) (2.1) 

9=0 

can be calculated iteratively in the instanton number, 1/a 4 powers below, using local- 
isation in the moduli space of gauge theory jT^l EH H3 HZ] or worldsheet instantons 
[HIE]- Global properties under the monodromy group, which should be central in the 
solution of the theory, remain obscure in these approaches. We find that studing the 
global properties enables us to solve the B-model completly and in particular fix the 
holomorphic ambiguity. This yields closed modular expressions which determine the 
instanton contributions to all orders in the instanton number, but are iteratively in 
the genus i.e. in A. We finish the section with some speculations how to obtain closed 
expressions in A as well. 



2.1 Modular properties of the genus zero and genus one sec- 
tor 

We focus on simplest case of SU(2) gauge theory without matter. Generalizations to 
other gauge groups and matter spectra with asymptotic freedom are certainly possible. 
The monodromy group of pure SU(2) Seiberg-Witten theory is T(2) G T = SL(2, Z) 
generated by [TH] 

Moo = PT- 2 , Mi = ST 2 S , (2.2) 
where S = ( J "J ) and T = ( J \ ) are the generators of F„ and P = 



Using modular properties as well as the B-model holomorphic anomaly 



-1 
-1 

we are able to express the J 7 ^ in terms of T(2) forms and the quasi- modular form of 
degree two E 2 . 

The natural embedding of gauge theory in type IIB string theory [20] is the ex- 
planation that the higher genus worldsheet technique of [I] applies to the space-time 
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instanton calculation. More precisely as found in [20] IIB theory on the local Calabi- 
Yau space 0(—2, —2) — > P 1 x P 1 has a double scaling limit in the two complexified 
Kahler parameters t±, t 2 of P 1 x P 1 in which !F^ 9 \ti, t 2 ) approaches F^'la). What we 
find below is that the holomorphic anomaly equation make sense in the limit and can 
be directly viewed as property of the gauge theory. 

The Picard-Fuchs equation for the periods a = J a X and od = A of the elliptic 
curve [1^] 

y 2 = (x-u)(x-A 2 )(x + A 2 ) (2.3) 
with meromorphic differential A = ^^—rdx 



2tt x 2 -1 

d 2 a 



d 2 u 4(1 - v? 



(2.4) 



o WO) 

allows to calculate the genus zero prepotential using the relation = a D up an 
irrelevant constant. Here for convenience we have set the Seiberg-Witten scale A = 1, 
which can be easily recovered by dimensional analysis. The details of the calculation 
that leads to explicit expressions for the periods can be found 1 in [21]. The elliptic 
curve (|2.Hjl with T(2) monodromy has the j-function j(r) = ^it^fl-i) 2 ' which allows 
to write [22] 

u{r) = , (2.5) 

as the Hauptmodul of T(2) in terms of the ratio r = f b —/ J a — — —4^^^ of the 
periods over the holomorphic one-form. The latter are solution of the 2 -Pi hyperge- 
ometric differential equation |^ = (jz^yf^ + inz^j 1 which implies that the inverse 
relation to (|2.5|) can be written in terms of a Schwarz-triangle function (see e.g. [2*3] ) 

We defined in (|2.5|) b = 9%, c = 9$, d = 9\ with the relation b + d = c. Using the 
modular transfomation 

h h b —> —r 2 d 

T: 7 S: d ^-r 2 b (2.7) 

C 4 * & 



C — > —T 2 C 



and (|2.2|) it is immediate that u(r) is invariant under T(2). For further reference let 
us note that the discriminant of the elliptic curve is given by u = 00 and 

A = u 2 - 1 = 4^ = . (2.8) 
b z 



1 In ~""I] the isogenous curve y 2 = (x 2 — it) 2 — A 4 with the meromorphic differential A = 1 ^2x 2 ^- 
is used. This means that a sw = a KLT , af^ = a^ iT /2, t^ w = t^ lt /2 and the monodromy is 
T (4) instead of T (2). 
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The Kahler potential in rigid special geometry is given by (see for a recent 
review) 

K = i^Fj - X J F/) . (2.9) 
For the Seiberg-Witten curve X = a and r := J^F so that the metric becomes 

g a - a = d a B- a K = 2r 2 . (2.10) 

The genus one amplitude is obtained by integrating the genus holomorphic anomaly 
using the special geometry relation 



=-log(Vw?) • (2.11) 

Note that the holomorphic ambiguity is fixed by requiring that F' 1 ' is invariant under 
SL(2, Z) transformations and regular inside the fundamental domain 7i. Indeed i] is 
the unique modular form of weight ~ so that regular inside Ti. 
In the holomorphic limit 2 f — ► oo, we get 

f« = _ log(r?(r)) = _ i log( ^) _ _L i og(M 2 _ 1} (2 12) 

This agrees with The form of JF*- 1 -' in terms of u and |^ follows in the rigid limit 
from and was observed in jS]. Using (|2.12j) and 77 12 = 2~~ 4 bcd one gets 

da 60 , 

m = t ■ < 2 - 13 ' 

Further with (j2.5|) and the formulas for the derivatives of the # functions we find 
-(o)/ s Or 8Vb 1 fdu x 3 



^'^^d^U 1 ■ (2 - 14) 



This can alternatively derived by taking twice the derivative w.r.t. to a in the Matone 

9a 



relation F (0) - §a^^ + w = [29 . Using the Picard-Fuchs (Q equation and (|2~T3j) 
we can write the period as' 
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a(r) = ^{E 2 + c + d) . (2.15) 



Closely related expressions for the SW-periods of elliptic curves appear in 14 and 
for more general gauge groups in [2*7] . 

The strong coupling duality element S : t ^ — ^ of iV = 4 Super- Yang-Mills 
theory is not trivially realized in pure N = 2 Seiberg-Witten theory. It relates 

2 Because of the volume of the diffeomorphism group generated by the Killing vector field on 
T 2 we need regularization of an infinite constant in that limit. Only the immediately well-defined 
d T F^ plays a role in the following. To define that limit consider r and f as independent variables. 

3 To keep the notation simpler we supress normalisation factors of ^7 from the periods, which 
make them an integral representation of (|2.2(1 . 
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the gauge instanton expansion in the region of asymptotic freedom -,- — > to 
the magnetic £7(1) which is weakly coupled at the magnetic monopole point z = 
(u — 1)/2,od — ► 0. The gauge coupling of the asymptotic free theory is determined 
by r(a), while the one of the magnetic U(l) is determined by td(cld) with the relation 

t d = — (2.16) 

T 

We note that the one- loop amplitude (|2.11|) is S'-duality invariant. The holomor- 
phic limit at the monopole points is taken by td — > oo. Therefore one has 

F D V>{a D ) = -logfafo)) = -^log(^) - ^log(« 2 - 1) 

log(2) i log(ajj) a D 3d 2 D 19a 3 D 4 (y2 ' l7 " > 

= 5 T7^ H U (an) 

3 6 12 2 5 2 9 2 12 3 V 7 

We derive, similar as ()2.15|) was obtained, from the first line of ()2.17|) a formula for 
the second period 



a D (r D ) = -^(E 2 -b-c) . (2.18) 



This can then be inverted to obtain the series expansion in the second line of ()2.17|) . 
where we rescale 

a D := i°f . (2.19) 

It is naturally to define an anholomorphic period 

1 

M 

with the anholomorphic weight one form 



Mr, f) = ^(E 2 {r, f) + c + d), (2.20) 



E 2 (t,T) :=E 2 (t) + — . (2.21) 

n(T — T) 

Different then the Eisenstein series E 4 and Eq, holomorphic modular forms of weight 
4 and 6 which generate the ring of holomorphic forms of T , the holomorphic E 2 
transforms not quite as modular form of weight 1, but rather quasi- modular under 
r , i.e. with a shift 

ar + b\ , . rt2 „ , , . 12 



E 2 — = (cr + dYE 2 (r) + —c(cr + d) . (2.22) 

\ct + a J 2,ixi 

However the anholomorphic piece in ()2.20j) cancels the shift so that E 2 (r,f) trans- 
forms as a form of weight 1. 

It follows from ()2.20|) that a(r) = lim^oo A(t, f) and a D {j D ) = lim^^oo ^A(—^, — -^). 
Morover S-duality transformations (|2.7|) allow us to write a and do as functions of r 
or Td- In particular we can integrate ()2.15j) w.r.t. ap in the magnetic phase 

?>«>> = l -al log (-f ) + a D - + f + ^2 + H# + 0(a%) . (2.23) 
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2.2 Propagators and integrating the holomorphic anomaly 
equations 

From the genus zero and genus one amplitude one can derive the propagator PQ. With 
the simplification for local B-model explained in [26J we obtain 



2F? If 3 \ 1 - 

sm = s =7r=4 B *M-^)= : 24 E ^>- (2 - 24) 

which in the holomorphic limit becomes 

2^ 1 

S aa = S= (2.25) 

This quantity contains the contribution from the boundaries of the WS moduli space 
in the topological string theory. It is closely related to Tj^ the quantity that arises 
if one considers correlators of integrated two-form operators constructed from the 
descent equation in the gauge theory on four manifolds. More precisely intersections 
of the correponding 2-cycles require contact terms which are Tjy = S aiCLj f^f^ 1 , see 



Equations (|2.14|2.1ip2.24jl define the data needed to recursively solve the 5-model 
PP. Using the fact that the formalism of integrating the holomorphic anomaly equa- 
tions commutes with the double scaling limit taken to obatin the gauge theory [2H] 
and power counting in the propagator in the Feynmann rules of pQ we obtain the 
following general result 

3(9-1) 

i^)(r,f) = X*- 1 4*(T,f)c?>(T) • (2-26) 

k=0 

Here we defined X = 172 g c 2 d 2 = io8(u 2 -i) 2 6 3 ' wn i cn transforms as a weight —3 object 
X(^) -> ^I(t) under T(2). is invariant under T(2), which implies that 

c£ are homogeneous of weight 3(g — 1) — k in (6, c, d). Further conditions on cjj? 
come from regularity at u = and a gap condition at the conifold u = 1 as dicussed 
below. 

We obtain the holomorphic limits of the expansion in the asymptotic freedom and 
the strong coupling region as 

T {9 \a) = lim F {9 \r,f), and T D{9 \a D ) = lim F {9 \- — , -— ) (2.27) 

and the a or an (&d) expansion are obtained by inverting (|2.15|) or (|2.18|) . We note 
that the leading behaviour of electric and magnetic expansion in these parameters is 

*"v~m^wF* and (2 - 28) 
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respectively. The first asymptotic behaviour can be derived in the gauge theory 
limit of type II theory on 0(— 2, — 2) — > P 1 x P 1 . More precisley one uses in the 
Gopakumar-Vafa expansion - — : 1 mX a = ^ Q A 2g ~ 2 (— l) (g-1) 2 ^ 3 _ 2 \\ m2g ~ 2 an d the 

1 2 sin — J y\ y )• 

multiplicity n~J = <5 9 ,o<5m,o = — 2 of BPS states corresponding to constant maps on 
one P 1 as well as properties of the limit discussed in [20] • The derivation is similar as 
for the constant map contribution Jj^-Cg_ l = 2g (2g-2)(2g-2y. m ^8J- ^ ne as y m Pt°ti c m 
the magnetic expansion comes from the occurrence of the c = 1 string at the conifold 

M- 

We come now to the explicite iterative solutions in the genus. For example the 
recursive definition of is 

F<»(t, f) = \ SF f + \s{F?f + L&iipf f - \lPFf - \^Ff + Xc%\u) , 

(2.29) 

were Cq (u) is the holomorphic ambiguity at genus two. This ambiguity must be 
invariant under T(2), which implies that it can be written in terms of u. Moreover 
ree; ularity of at u — > oo and the leading pole behaviour at u — > 1 implies that it 
is of the form 

2 9- 2 A (g) 

X°-^\u) = u*-°J2^ (2.30) 

i=i ^ > 

Note that A\ are undetermined constants and the right hand is a rational function 
of the r(2) invariant function u. Using (j2.14|2.11j) and the standard formulas for the 
derivatives of 9i, E 2 we may write ()2.29j) with h = (b + 2d) as 

F>»(r, f) = f (SB| - 9Elh + 6E 2 (f + cd) - Ml^lliM ) , (2 . 31) 

an almost holomorphic modular function of T(2). Here we determined the ambigu- 
ity as follows. Using (j2.27l2.15j) we can expand ()2.31j) in the electric and magnetic 
holomorphic limits. With the leading behaviour ()2.28|) we found A 2 = — 12 1 9 9 60 and 
A {2) = - Jg. This yields 

^ = - ^ - - ^ + O(i) (2,2) 



240 -2 a 2 2 18 a 10 2 22 a 14 2 34 a 18 V 2 ' 

This series predicts all genus 2 instantons and checks with the coefficients that ap- 
pear in the literature JH]!B]- The expansion in ao is obtained from (j2.31j) using 
(I2.27l2.7l2.l8l) 

F D{2) = K - % - ^ - + 0(5 4 D ) (2.33) 

240 a 2 D 2 13 2 16 2 17 5 V d; V ; 
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Solving the recursion for genus 3 yields 

F ( 3) = SF W F W _ isp F M F f + l SF f + \S\F^fFf - \S*Ff{F { ?Y 
-\S\F { l ] f{Fff + '\S\F^fFf + SlF^F^Ff - i&F^Ff 
-\S\F^f + \S 5 Ff{Fff - IS'F^FfFf - \s^Ff{Fff 

i 1 C3 7^(1) 7^(0) 1 5 Q3 p(0) , 1 r-3 77(0) pft) ' 1 r-2 ' 7 r-4 p(0) p (0) 

+ 4^ ^,2 ^,4 +12^ ^,3 ^,3 +8^ ^,5 ^,1 "8^ b A ~ 48^ ^,5 ^,3 

+ I^^4 0) (^.3 0) ) 2 - ^^(4 0) ) 4 - ^ 4 (Ff ) 2 + h^ F f + ^\u) , 

(2.34) 

and we determined the coefficients = 244 5 9 9 440 , A^ = g^fffo, = 25ifio anc ^ 
= 3§|f|o- This yields the following almost complex modular expression for F^ 

F (3) = 4X 2 ^ 5Ej 6 _ 2hElh + 40£ 4 (26 2 + hcd) - §£ 3 /i(529& 2 + 559crf) + 

(11726 4 + A0Q0b 2 cd + 1223(cd) 2 ) - ^f/i(8446 4 + lQ85b 2 cd + 310(crf) 2 ) + 

2Yq(107186 6 + 495966 4 cd + 440076 2 (cd) 2 + 944(cd) 3 ) 



(2.35) 



from which the electric 

1 ( \ 441 18459 62106849 256368735 1 . 

= 1 1 1 1 h 0( — 

1008 -2 3 V a 4 2 12 a 12 2 16 a 16 2 28 a 20 2 31 a 24 l a 28 



(2.36) 



and magnetic expansions 



D{3) _ . . 9a D 143 638275^ 
^ "T0085 4 ; "2^~ 227? + (A^J (2.30 

follow from ()2.27|) . The number of terms in the modular expressions of F^ grow much 
slower then the number of graphs in the holomorphic anomaly equation, because many 
graph contributions are proportional to the same quasimodular form. For genus four, 
where the holomorphic anomaly equation has 83 graphs we find 

F (4) = 4X 3 {^-E^- 985Elh + 2E 7 2 (2399b 2 + 7001cd) 

-— E%h(3761b 2 + 6125cd) + -£ 2 5 (858636 4 + 3633446 2 cd + 240083(crf) 2 ) 
3 5 

--E 4 /i(6044696 4 + 16773406 2 cd + 547811(cd) 2 ) 
5 

+ ^£ 2 3 (5312666 6 + 27936156 4 crf + 32851236 2 c 2 d 2 + 447656(cd) 3 )) 

— E 2 /i(44307566 6 + 165503376 4 crf + 119259276 2 (cd) 2 + 889964(cd) 3 ) 
35 

1 - 

+ £ 2 (312324286 8 + 1952748406 6 crf + 3296138196 4 (cd) 2 + 130729960& 2 (cd) 3 

L \ 

+3566728 (cd) 4 ) - y^r/i(87826748& 8 + 4237709486 6 crf 

+5113136016 4 (cc/) 2 + 1280981 726 2 (cd) 3 + 4442006 (erf) 4 )} (2.38) 
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with electric 

, A , - * ( 1 765 126195 1925006715 14420664765 m , 1 , 

tW — —1440 .9 I I I I I- O 

1 I a 6 2 12 a 14 2 16 a 18 2 29 a 22 2 32 a 26 l a 30 ' 



(2.39) 



and magnetic expansion 



m 1 1125 ~a D 3915 a 2 , 4786021 o| 

^ ~ ~1440^ + ^ 2^ + 2^ +U{aD) {2AU) 

We derived expressions for in terms of modular forms up to genus six and 
checked the large a expansions against results made available to us by Nakajima 4 . 
Let us report here the dual expansions, which are interesting as they correspond to 
perturbations of the c = 1 string at the selfdual radius by momentum operators 

t d® _ 1 77175 ~a D 100971 ^ 51425582135^ . 4 ( . 

* " 10564 "^36 2^ 2«U + U ^> ^ 



D(6) _ 691 187535255 D 169085255^ _ 672990085791^ 4 

~ "327600 5f 2 s 2^ 24913 + ™ ' 

(2.42) 

2.3 Fixing the ambiguity 

In ()2.26|) all cjj. except for Cg , the holomorphic ambiguity, are determined by the 
recursion relations ()2.29I2.34|) that follow from the holomorphic anomaly equation 
in terms of lower genus F^ g \ Genrally fixing the holomorphic ambiguity is a major 
problem in the B-model, which however the case at hand is completly solvable. The 
discriminant locus of the curve ()2.3j) is at u = oo and at u — ±1 where ()2.3)1 develops 
a node. At all other points in the moduli space of the Seiberg-Witten geometry F^) 
must be regular. As follows from the global properties of the 9 functions and E2, 
regularity at u = restricts the form of cf^ (r) to 



c 



(9) 

k 



k 

(l + (_l)9 + fc) 



h s P d(9 , k) (b\cd), (2.43) 



where Pd( g ,k)(b 2 , cd) is an homogeneous polynomial in b, c, d of degree 

d(g, k) = 3(g - 1) - k - { \ } ) , (2.44) 



with d ^' fc - > + 1 coefficients. In particular in the ambiguity the number of unknown 
coefficients of this polynomial £ 
(I2.30|) . which grows with ~ 2g. 



coefficients of this polynomial grows with ~ |<? slower then the number of the Af in 



1 These somewhat lengthy expressions are available on request. 
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Moreover the leading terms in (|2. 3312. 3712. 40p2.4ip2.42JI correspond to correlators 
of the cosmological constant operator of the c = 1 matrix model in the genus g vac- 
uum sector j3U]. Very important is the occurence of the gap in J r ^ D , i.e. the absence 
of terms k = 0, . . . , 2g — 3. The 2g — 2 gap conditions fix the Af , i — 1, . . . , 2g — 2 

a D 

constants and hence the ambiguity ()2.30|) . Together with the anomaly equation this 
provides a very efficient way to solve the theory completly. The asymptotic ()2.28j) 
and the particular form of 1)2.43)1 are further consistency constraints confirming the 
gap property. A gap follows if there is a matrix model expansion for the holomor- 
phic topological string at a critical point, as e.g. at the orbifold point in the local 
0(— 2, —2) — > P 1 x P 1 model discussed [7]. In this case the measure integration yields 
at each genus the negative power term in expansion parameter and the perturbative 
terms start with positive powers. This particular behaviour of the magnetic expansion 
of the N=2 pure SU(2) model at the conifold is hence explainable by the proposed 
unitary matrix model |2]. As we saw above the model is solved by the gap property 
and the holomorphic anomaly. If one would like to employ techniques in order to 
determine the weak coupling instanton expansion one would have to do exactly what 
we have done in ()2.26j) namely to write the result globally. More generally the gap 
can be understood from the absence of correlators of the ground ring operators [20] in 
the c = 1 string describing the limit of the toplogical string near the nodal singularity 
(conifolds) of local models. Indeed we have checked that the gap occurs also at the 
conifold in the local 0{— 3) — > P 2 geometry and fixes the ambiguity of this model. 
How this extends generally to singularities of local models and the modifications for 
singularities of global models will be discussed in |5T] . 

If we absorb X into genus expansion parameter A in (J2.26J) . it becomes a sum of 
a quasi-modular form. The simplest example of such an expansion, where the coef- 
ficients are however modular forms of the full modular group Tq appears in Hurwitz 
theory on T 2 [32] • As reviewed there this leads directly to combinatorial problem 
that is solved by free fermions and Z = e F can be written in a product form that has 
been recognized as a generalized ^-function product form in [22]- More examples of 
such product forms are provided by vertex algebras [211 EE] and arise in heterotic type 
II duality [2U [23 I2H1 123- It would be interesting to see whether the SW partition 
function is related to a generalized T(2) theta function. 

3 Dijkgraaf-Vafa conjecture 

In jH], Dijkgraaf and Vafa proposed a remarkable relation between B-model topological 
string on a non-compact Calabi-Yau geometry and a matrix model. 
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3.1 Dijkgraaf-Vafa transition and geometric engineering 

The n cut matrix model is obtained by reducing holomorphic Chern-Simons theory 
on D5-branes wrapping n P^s in a modification of the geometry 0(— 2) © O(0) — > P . 
The k-th P 1 is wrapped by branes, k = 1, . . . ,n. In the modified geometry the 
location of the P^s in the originally flat O(0) x-direction is now fixed at the minima 
of a potential W(x) of degree n + 1. E.g. the n — 1 geometry is the blown up conifold 
O(-l) © O(-l) — ► P 1 . The reduction yields a complex bosonic matrix model with 
the matrix potential 5 W'(x), that needs as additional data the choice of contour for 
the eigenvalue integration 0. 

The B- model geometry emerges after a transitions in which the n P x 's are shrunken 
and deformed to S 3, s. It has a local description as a hypersurface in C 4 

vw = W'(xf + f(x) + y 2 (3.45) 

where x,y,v,w are coordinates of C 4 , f(x) is polynomial of degree n — 1 that splits 
the n double zeros of W'(x) 2 . 

The latter geometry has been considered in to geometrically engineer Af = 1 
four- dimensional supersymmetric gauge theory. After the transition the breaking to 
M = 1 is achieved by putting units of Ramond flux on the S 3, s and the topological 
string or the matrix model calculates terms in the Af = 1 effective gauge theory. In 
|Hj it is already shown that the special geometry relation that determine the tree 
level (genus zero) topological string amplitude on the geometry ()3.45|) arise from 
the planar diagrams of the corresponding matrix model. The planar loop equation 
of the matrix model gives the spectral curve of the local geometry and the effective 
superpotential in a large class of Af = 1 gauge theory can be computed exactly by 
the genus zero amplitude of matrix model or topological string. It is conjectured 
that higher genus topological B-model string amplitudes should also be computed by 
higher genus diagrams in the matrix model. 

The meaning of the topological string amplitudes F( 9>0 \Si) in the effective the- 
ory is as follows. In Af = 2 supergravity action they determine the exact moduli 
dependence of the F-terms 

J d A x J d A eF^(Si)(W aP W aP ) 9 (3.46) 

here W Q/ 3 is the M = 2 graviphoton superfield and the Si are the M = 2 vector 
multiplets whose complex scalar field corresponds to the moduli of Calabi-Yau man- 
ifold. After integrating over the Af = 2 superspace these terms become the coupling 
-R+-F+ s ~ 2 of the self-dual part of the Ricci tensor R + to the selfdual part of the 
graviphoton field strength F + . 

5 Further generalization of this conjecture to the case of Calabi-Yau geometry with ADE type 
singularities can also be made 
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After breaking of J\f = 2 to Af = 1 by fluxes the topological string amplitudes 
F^ 9>0 \Si) occur in the following two terms of the Af = 1 action f4Tl 1321 H3j 

r\ = gj d 4 xj cPeWa^W^iFstF^y^F^iSi), (3.47) 
d'x J d 2 9 (F a pF a P) 9 Ni ~~^r~^ • (3-48) 



Here W Q/ 3 7 are now Af = 1 gravitino multiplet, and Sf are the Af = 1 glueball chiral 
superfields coming from the original Af = 2 vector multiplets jH]. The graviphoton 
field F a p can be treated as background field in the Af = 1 theory. In [UJ a C- 
deformation is introduced to deform the anti-commutation relation of gluino ip a to 
the followings 

{ij a ,^} = 2F af3 (3.49) 

It is shown that the effect of turning on the graviphoton background F a p can be 
captured by this C-deformation, and the F^ in the second contribution T 2 (|3.48J) 
are computed by matrix models at genus g. We can also see that the first term Ti 
in ()3.47|) contributes at genus g = 1 even there is no graviphoton background F a p. It 
is shown in [12] that this genus one contribution is also computed by matrix model 
genus one diagrams. There are also some other types of gravitational corrections 
besides (j3.47p3.48)) of the form Wa^W^ 1 S n from planar diagrams, which become 
trivial after the extremization of the glue ball superfield S jlH]- Our results confirm 
these very interesting ideas in [3TJ 021 HH] by a first direct tests of the connection 
between topological strings and matrix models at higher genus without using the 
superspace techniques in the effective gauge theory. 

3.2 The two cut geometry and the tree level and genus one 
amplitudes 

We consider now the case of a cubic potential W(x) = |i 2 + |x 3 in the Dijkgraaf- 
Vafa geometry (j3.45j) . The degree one polynomial f(x) = fiix + /io splits the double 
zeros of W (x) 2 at x — ai and x = a 2 to the four roots af, af of the equation 

W (x) 2 + f(x) = (3.50) 

We adopt the notation of [HI El, and change variable (a±, af, aj, af) = (xi, x 2 , x 3 , x 4 ) — ; 

(zi,z 2 ,Q, I) where 

Zl = ^(x 2 - xi) 2 , z 2 = ^(x 4 - x 3 ) 2 

„ 1 , .TO 

Q = -zixx +x 2 + x 3 + x 4 ) = 

^ 9 

I = I [( X3 + X4 ) _ ( Xl + X2 )]2 = -2{z x - z 2 ) . (3.51) 
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The local Calabi-Yau geometry ()3.50|) depends really only two complex structure 
deformations. We will use (zx,z 2 ) or below the A-periods Si and S 2 to parameterize 
them. The dependence of the genus zero super- and the higher genus potentials 
can be reconstructed from scaling laws and dimensional considerations. In particular 
we frequently set g,m = 1. 

The fundamental periods of the local geometry are 

i r a t i f A 

Si=o- A ' n i = — / A , 
2vr2 J a - 2m J a + 

where i = 1,2 and A = dx^JW 1 '(x) 2 + f{x) is a meromorphic differential which 
emerges after integrating the holomorphic Calabi-Yau (3,0)-form over an S 2 fibre 
direction of the S 3 's in (I3.45j) [3] . In [I] the integrals where further calculated pertur- 
batively for small Z{. This limit corresponds to vanishing S 3 s and is suitable for the 
perturbative matrix model expansion. Solving the B-model and fixing its ambiguity 
requires a global understanding of the complex moduli space in z 2 ). We therefore 
derive the Picard-Fuchs equations and use them to explore the global properties of 
the integrals Si and IT. 

We find that derivatives of A w.r.t. Zj up to second order multiplied with suitable 
polynomials in z% combine to a total derivative, i.e. £«A = J^k+^P^'^i^)^^^ = 
dgi(x, zi, z 2 ). Naively the differential ideal with this property is generated by three 
independent differential operators £j. However dgi(x, Z\, z 2 ) is a meromorphic differ- 
ential with non- vanishing residua, hence one cannot conclude from the exactness that 
d J r A = 0. For the following two operators the residua vanish 

£i = 2zi(2*i + 6z 2 - l)d 2 Zl + (1 - 10zi + !2 Zl 2 + A Zl z 2 )d zl d Z2 + (3 - 2z x - 6z 2 )d Zl 
+ (1 <-> 2) 

C 2 = (2 Zl + 2z 2 - l)[-2zi(l - A( Zl + z 2 ) + W - 2z x z 2 - W)d\ x + {z x + z 2 ) 

(1 - 8 Zl + Q Zl 2 - Qz lZ2 )d Zl d Z2 ] + ( Zl (7 - 18z 2 + 2Qz 2 2 + AQ Zl 2 + Zl (Q2z 2 - 39)) 
-3z 2 (l - 3z 2 + 2z 2 2 )))d Zl - 3(1 - \2 Zl + 18^ 2 + U Zl z 2 ) + (1^2). 

(3.52) 

These Picard-Fuchs operators annihilate the periods and fix their expansion up to 
linear combinations. The discriminant of these differential operators has the following 
components 

Ci : Zl = 0, C 2 : z 2 = 0, I : / = 1 - 2z x - 2z 2 = 0, 
J: J = (xi — x 2 ){x 2 — x 3 )(x 2 — x 4 )(xi — x 4 ) (3.53) 
= (1 - 3 Zl ) 2 - Qz 2 + 9z 2 + Uz x z 2 = 0, 

whose schematic intersection after a suitable desingularisation of three order tangen- 
cies is depicted in Fig. 1. 
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Figure 1: Divisors in the moduli space of B-model, (zi,z 2 ) = (0,0) is the matrix 
model expansion point and along the dash line one has enhanced N = 2 SUSY. 

According to the Dijkgraaf-Vafa correspondence the periods Si are identified with 
the filling fractions Ni, with Y^i=i ^ = ^, of eigenvalues in the large N limit 
of the dual matrix model, and it is shown that the special geometry relation and 
Picard-Fuchs equations are reproduced in the planar limit of the matrix model [3] 
and the genus zero topological string amplitude follows from integrating the 
special geometry relation = n*. To analyze the exact effective superpotential 
W e ff = 27ri^ i (A^ i n i + ctjSi), where N and cti are 3-form flux quanta through the 
Ai and Bi cycles respectively, globally we need to the periods troughout the complex 
moduli space. This is done in appendix A, where we also find that there is no point 
where the periods degenerate quadratically in the logarithms. Which is the signal of 
a large volume or maximal unipotent point in the moduli space. 

Also at one loop the conjecture holds 6 (HI CD]. The B-model expression for the one 
loop free energy obtained by integrating (|4.59jl using ()4.58|) and fixing integration 
constants it turns out to be [Hj 

F {1) = \ log (det (J|) (z lZ2 )-hj^ , (3.54) 

where we obtain the periods Si in terms of complex structure moduli Zi and the 
corresponding inverse series 

z 2 ,g) = \z x - ^z 1 (2z 1 + 3z 2 ) - f-Zi(4z 2 + 13 Zl z 2 + 9zj) + 0(g 6 ) 
S 2 (z 1 ,z 2 ,g) = -S l {z 2 ,z 1 ,g) 

Zx{S x , S 2 , g) = ASx + 8g 2 S 1 (2S 1 - 3S 2 ) + 8# 4 (20S 2 - 67S& + 39S 2 2 ) + 0(g 6 ) 
z 2 {Si,S 2 ,g) = z 1 (-S 2 ,—S 1 ,g) (3.55) 

6 The one loop-test in^Hl and the higher loops tests [Hj are made at the N — 2 point Si = — S^- 
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from the two power series solutions to ()3.52j) at (zi,z 2 ) = (0,0). Identifying 5, 
with the filling fractions iVj we get the genus one contributions to (|B. 101(1 . Note 
that the integration constants c$ in (ziZ2) Cl F 2 J° 3 , which fix the behaviour of 
at the discriminant components are global data do not depend on the base point 
(zi,Z2) = (0,0) or the holomorphic limit Si — > taken at this base point to obtain 
the matrix model expansion. The coefficient c\ — — | at the conifold (shrinking S 3 ) 
is a universal property of the topological B-model. 

To solve the B-model recursion we need the genus zero three point functions, 
which are rational functions in complex structure moduli. 

1 - (6zi + 5z 2 )g 2 + 3{3z 2 + 3z x z 2 + 2z 2 )g 4 1 - {3z L + 5z 2 )g 2 

212121 lQg^P ' 212122 lQg 2 I 2 

(3.56) 

The three point functions C ZiZjZk are symmetric in ijk and from symmetry consider- 
ation follows C Z2Z2Z2 = C zxzxzx {z x <-> z 2 ) as well as C Z2Z1Zl = C zlZlZl {z x «-> z 2 ). 

The corresponding matrix model is a Hermitian matrix model with the cubic 
potential W^($) = |$ 2 + |$ 3 for a rank N Hermitian matrix $. The partition 
function and free energy F of the model are 

In the large N limit the eigenvalues distribute around the two critical points and 
— i of the potential and form two cuts. We consider the metastable vacuum where 
with N\ eigenvalues at and ./V2 eigenvalues at — -. This is a two-cut solution of the 
matrix model with N\ and N 2 fixed and subject to the condition N\ + N 2 = N. In the 
large iV limit the free energy of the matrix model has genus expansion in 1/N 2 , and 
at each genus there is a perturbative expansion by the t'Hooft coupling constant gN. 
Dijkgraaf and Vafa conjecture that the free energy of the matrix model at each genus 
is matched to the topological string amplitudes on the local Calabi-Yau geometry 
(|3.45|) . by identifying of the periods Si in the geometry with the eigenvalue filling 
fractions N{. In Appendix El we review more details of the matrix model calculations 
of the free energy. 



4 The holomorphic anomaly equations 

The key to the solution of the topological B-model are the holomorphic anomaly 
equations. To solve them recursively one needs in general to derive three types of 
propagators S t: > , S 1 ^ and S^ [Tj. In local geometries S 1 ^ and S^ can be gauged to 
zero ^26] and the derivation of the propagators in the multimoduli case is discussed 

in mcsiEsj. 
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4.1 The holomorphic anomaly recursions 

The geometry on the complex structure moduli space Zi of Calabi-Yau is a spe- 
cial Kahler geometry. Its metric, connection and curvature are determined by the 
Kahler potential K by the well-known formula, Gq = didjK, Y\ m = —G lk diG km and 
R k -., = — djTft. They have a well-known special geometry relation with the three point 
Yukawa coupling Cijk = DiDjD k F^°\ which comes from the tt* equation and can 
be thought of as the holomorphic anomaly equation at genus zero 

4* = + Gtfi ~ CaJCf" (4.58) 

At genus one and higher genus the topological string amplitudes has a holomorphic 
anomaly, and the anti-holomorphic dependence of the genus g free energy is related 
to lower genus free energy by the holomorphic anomaly equation pQ 

didjFU = \c m C kl - - i) Gil 

1 _ 9-i (4.59) 
diFb) = -CtiDjDkF^ + J2 D i Fir)D kF {9 ~ r) ) g > 2 . 

2 ' r=l 

Here the Cf = e 2K C- tfk G^G kk define the propagators S ij as diS jk = G k . The 
holomorphic equation can be integrated and represented as graphic Feynman rules 
to give the higher genus free energy in terms of lower genus up to a holomorphic 
ambiguity. The propagators can be solved by integrating its defining relation and use 
the special geometry relation (|4.58j) . One finds 

S*C m = 5\d k K + + r kl + f kl (4.60) 

here f % kl are ambiguous integration constants, and they are meromorphic rational 
functions of the complex structure moduli Zi with poles at discriminant points of 
the moduli space. Suppose there are n complex structure moduli, then there are 
|n 2 (n + 1) equations for \n{n + 1) propagators S l K In the case of one modulus, the 
number of equations and propagators are the same, so the meromorphic functions /j fc 
can be just set to zero. However, in the multi-moduli case we consider, the equations 
over determine the propagators, so we have to choose the ambiguity /j fe properly to 
satisfy some constrains and ensure we can solve for the propagators. 

There are certain simplification in B-model calculations for the case of non- 
compact local Calabi-Yau manifold. In this case there is a choice of gauge such that 
the Kahler potential K and metric over the moduli space Gg in the Si coordinates 
is a constant in holomorphic limit, and the dilaton component in the propagators 
vanish. So in the holomorphic limit Sj — > the connection vanishes and covari- 
ant derivative in the Si coordinates is just ordinary derivative. This also makes the 
topological string amplitudes entirely independent of quantities such as Euler number, 
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Chern classes of the Calabi-Yau, which will need to be regularized in the non-compact 
case. In this case, the metric and the connection in the z\ coordinates are 

dS k 



dz, L 
dz. 



C, 



d 2 S, 



ZjZ k 



dSi dzjdzk 



(4.61) 



where C k j are constant in the holomorphic limit. 



4.2 Propagators and Dijkgraaf-Vafa conjecture at higher genus 

The geometry we consider has two complex structure moduli. In order to solve the 
propagators, the holomorphic ambiguity /j fc have to satisfy 3 constrain equations 
by eliminating propagators in (|4.6U|) . These constrains equations are rational 
functions of the complex structure moduli Zi, 7 and we are able to find a rational 
solution for the 



fh 

fl2 

fh. 



6 - (49zi + A8z 2 )g 2 + (163z? + 219z l z 2 + I26z^)g 4 
+(210z? + 30Az 2 z 2 + 242z x z 2 + 108^ )# 6 )] /(20g 6 z 1 I 2 J) 
- 29 - (79zi + I57z 2 )g 2 + (10z 2 + 260 Zl z 2 + 2l0z%)g i }/(20g 4 I 2 J) (4.62) 
7 - (552i + 68z 2 )g 2 + (142^ 2 + 315ziz 2 + 219^ 2 )# 4 
-(120z? + 338z 2 z 2 + 492zxz 2 + 234zf )# 6 ] / (20g 6 z 2 I 2 J) . 



By definition /L = f k j and the (zi <-> z 2 ) symmetry determines the other /- fc , e.g. 
/| 2 = fli{z\ <-» z 2 ) e.t.c. Note that the discriminant factors in ()3.53j) should be 
the only singularities appearing in the denominator for the ansatz of holomorphic 
ambiguities and in fact all appear. 

The holomorphic anomaly equation at genus one can be integrated to give the 
Ray-Singer torsion of the target manifold. The genus one free energy can also be 
expressed in terms of genus zero three point functions and the propopogators as 
follows 

d . F W = -S> k C l]k + ^ a r log( A r ) (4.63) 

r 

here a r are constants and A r are various discriminants of the local geometry. The 
holomorphic ambiguity /j fe should give a solution of the propagators S 1 ^ that satisfies 
the above consistency check ([4.63)1 . We have chosen our ansatz of the holomorphic 
ambiguity ()4.62j) that satisfies the (j4.63J) with the constants a\ = 1/15, a 2 = 2/15 



d t F^ = ~S jk C ijk + $(1 log(*i* 2 ) + ^ log(A 2 )) 



(4.64) 



7 We note that while genus zero three point functions are rational functions of the complex 
structure moduli 2$, the connections T l - k are generally not rational functions. They combine to 
give rise to rational equations for holomorphic ambiguity /j fe . 
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This choice of ansatz is convenient in the sense that it leads to the correct leading 
behavior in genus two, so it is easier for us to fix the genus two holomorphic ambiguity 
there. 

In local geometry the genus two topological free energy can be integrated from 
the holomorphic anomaly equation. It is 

F m = (4.65) 

+^S kl S™F^ m Ff n + l^S kl S mn F$F^ n + /CO 

We fix the genus two holomorphic ambiguity with some initial data from matrix 
model calculations 



/( 2 ) = - 1253 - 10503(zi + z 2 )# 2 + 27(10812? + 95(kiz 2 + 1081z?k 4 

L " (4.66) 

-26865(2! + z 2 )(zi - 2 2 ) 2 ( 7 6 ]/(9000(7 4 2 1 2 2 J 2 ), 

Once the holomorphic ambiguity is fixed, we can compute the genus two free 
energy to very high order using ([4.65)1 . The topological string approach is much more 
advantageous than direct matrix model calculations where it is hard to compute to 
free energy at higher orders (see Appendix El for more details). After some extensive 
computer running time, we are able to make many checks of the topological string 
predictions (|4.67|) below for the genus two free energy. 



F (2) = -— (-^ + -^) + (—N 1 -—N 2 )g« 
240 Vi 2 N 2 2> V 6 6 Jy 

66132 

+ (338iVi 2 - 1632AyV 2 + 338N 2 2 )g 8 + (^— iV\ 3 - 12088(Wi 2 iV 2 + • • • )g 10 



+ ,1305280^ _ 18059582^3^ + 11824166JVl2jV22 _ . . . )g u 

o o 

+ (12963696iVi 5 - 244438427iVi 4 iV 2 + 745362156iVi 3 iV 2 2 )g 

ATfi 35002369227, r - T 148144671957 , r . . 
+ (362264064iV 1 6 N^N 2 + N^N 2 



ii 



-57597548553iV 1 3 iV 2 3 + ■ • • )g w 

,29035470208 7 862430780350 fi , 9 

+ ( Ni iVi 6 iV 2 + 1580964252892iV 1 5 iV 2 2 

10256675032550 Ar4Ar o , 1R 
- N x N 2 +■■■ )g 18 

, 1250634104832 %r „ 44363662176978 , r - %r 303466060570354 %r , Ar 2 
+ ( Ni A^i 7 A^ 2 + Ni N 2 

Do O 

851102001082126^5^3 + 237 637137780236iV 1 4 iV 2 4 - ■ • -)g 20 (4.67) 



5 

The main difficulty in the B-model calculations is to fix the holomorphic ambi- 
guities at each genus. Also the Feynman rules that solve the holomorphic equations 
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quickly become very complicated. Here we push the calculations only to genus two 
in our calculations, since there are less new conceptual issues beyond that. 

5 Conclusion 

The B-model iteration in the genus bears some resemblance to the procedure compute 
higher genus free energy and resolvent of the matrix models for one-cut solution in 
[T2] and generalized to multi-cut solution in [121 EH1 EH] , where the iteration equation 
is obtained by doing 1/N expansion in the loop equations, and looks similar to the 
holomorphic anomaly equation in topological strings. 

From the hermitian matrix model point of view the anti-holomorphicity is very 
unnatural. The holomorphic anomaly equations were re-interpreted in [2] as infinites- 
imal manifestation of the fact that the topological string partition function transforms 
as a wave function under change of polarisation in the middle cohomology of the tar- 
get space. Using this picture the failure of holomorphicity can be traded against a 
failure of modularity with a similar iteration [SU], which makes the connection more 
naturally. 

It would be very interesting to compare this latter iteration with the iterations 
|48jl4*9] in detail, since this can in principle fix the holomorphic ambiguity in B-model 
calculations. Fixing the holomorphic ambiguity systematically is one of the main 
difficulties for topological string calculations in many other models. We hope further 
studies will clarify these issues and provide valuable lessons in fixing holomorphic 
ambiguity in more general models. 
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A Moduli space and monodromy of the two cut 
matrix model 

A.l Compactification of the moduli space and local expan- 
sions 

The aim of this section is to obtain the periods everywhere in the moduli space and 
to determine the monodromies. For the compactification of the moduli space we use 
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the projective space CP 2 with homogeneous coordinates (zi,z 2 ,z 3 ) and identify the 
z 3 7^ patch with 

zi = ^, z 2 = ^ . (A.68) 

Z3 z 3 

In addition to the divisors listed in (|H.5H|) we get now a CP 1 divisor at infinity, at 
which the periods turn out to be non-singular. We calculated the local expansion 
near all normal crossing divisors and determined the local monodromy. By analytic 
continuation we determined the global mondromy. One remarkable aspect of the 
geometry is that there is no point in the moduli space where at least one of the 
periods degenerates with double logarithm, which would correspond to the normal 
large complex structure point of a local geometry at which the mirror expansion in 
the large Kahler coordinates leads to a convergent instanton sum. 

Suppose we expand the Picard-Fuchs equation around a common point of two 
singular divisors Ai(zi, z 2 ) = and A 2 (zi, z 2 ) = 0. In order to find complete solutions 
of the Picard-Fuchs equation, one must choose a good local coordinate around these 
singular points. The technique for choosing good local coordinates is quite standard 
in algebraic geometry. For our two parameter model, there are two possible cases: 

• det(^- i ) 7^ 0, then the point A i (z 1 , z 2 ) = A 2 (zx,z 2 ) = is called the point of 
normal intersection of divisor Ai = and A 2 = 0. In this choice of good 
local coordinates is simply (Ai, A 2 ). 

• det(^ i ) = 0, then this is called a point of tangency of divisor A x = and 
A 2 = 0. In this case one will not be able to find all solutions around the point 
of tangency with the choice of local coordinates (A 1; A 2 ). We will encounter a 
very common situation in which the divisors have the following form 

Ai = a 2 + be 

A 2 = b (A.69) 

here a, b and c are degree one polynomial of complex structure moduli z\ and z 2 . 
The standard technique in algebraic geometry is to introduce two exceptional 
divisors to resolve the point of tangency. It turns out that a choice of good local 
coordinate in this case is (a, -%). In our analysis we will follow the standard 
procedure and use this good local coordinates. 

We list the asymptotic solutions of Picard-Fuchs equations and their monodromy 
at various singular points of the divisors (|3.5Hj) in the moduli space. Some of the 
singular points can be obtained by exchanging zi and z 2 , and we only list once these 
symmetric singular points. 

1. G\ flC 2 . his is the matrix model point we have tested Dijkgraaf-Vafa conjecture 
at higher genus. Th good choice of local coordinates is simply (zi,z 2 ). For 
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completeness we list the periods 

h = j -^1(2*1 + 3*2) + ••• 

1 z 1 

/ 3 = f 1 \og(z 1 ) + — + j-—(-4zj + z 1 z 2 + 5z 2 2 ) + --- 

U = f 2 \og(z 2 ) - ^ - ^ - ^(5zf + Zl z 2 - Az 2 2 ) + ■ ■ ■ (A.70) 

2. Ci fl J. The intersection point is at (-21,-22) = (0, |), and the choice of good 
local coordinates is (xi,x 2 ) = (zi, 1 — 2z 1 — 2z 2 ). The asymptotic solutions for 
periods are 

fi = y/x^il - 4xi - x 2 ) 

2 ^ "^2 1 



/ 2 = y/x 2 '(x 1 + — - (2xj - -X1X2 - -g-) H ) 

5 3 

/ 3 = xi(l--xi--x 2 + ---) 

/ 4 = / 3 log(xi) + ! + X2 + (^-3xix 2 -^) + --. (A.71) 

3. Ci fl J. This is a point of tangency of the two divisors at (z 1: z 2 ) = (0, |). The 
singular factor J can be written as 

J = 1 - Q Zl - 6z 2 + 9z\ + U Zl z 2 + 9z 2 2 = 9(z 2 - \) 2 + Zl (-6 + Uz 2 + 9 Zl ) (A.72) 

Following our discussion above we see a good choice coordinates is (xi,x 2 ) = 
{-, — ^tt^, k — z 2 ). This is the local coordinates around the intersection of the 
blow up divisor with the divisor C±. The asymptotic expansion for the periods 
are 

l/i 1 23 5 2 , X 

fi = Xix4{l — — — xi — — x\x 2 — 77^^X0 + ■■■) 

2K 216 72 46656 1 

4 

h = xl(l + -x 1 + 4x 2 + ---) 

f 3 = 1 + 81x 3 2 - 108xia;l - ™x\ + ■■■ 

U = /ilog(xi)+4(y-^2 + ---) (A.73) 

We can also solve the Picard-Fuchs equations with the local coordinates around 
the intersection of the blow up divisor with the singular divisor J. This will be 
useful later on when we try to match the basis and derive the monodromy of 
the divisor J. The good choice of local coordinates around this point is 

z l (Q-Uz 2 -9z l ) 1 
Xl = 9(z 2 - W ' X2 = 3" Z2 (A ' 74) 
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We find the asymptotic expansion for the periods with this coordinates 

j- 2 1/1 5 33 \ 

h = X\X% (1 - -r^xi + -^x 2 H ) 

lb 8 

2/1 3 43 N 

f 2 = x 2 [l + -x 1 - —X 2 + ---) 
h = 1-648x1 + .-. 

U = / 1 log(x 1 )+x|(^ + 32x 1 -^x 2 + ...) (A.75) 
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l i- 

write the singular factor J as 



4. I (1 J. This is a point of tangency between the divisors at (zi, z 2 ) = (|, |). We 



J = (zi - z 2 f + (1 - 2zi - 2z 2 )(l - Az x - Az 2 ) (A.76) 

l-2z 2 -2; 



A good choice of local coordinates is (xi,x 2 ) = ((zi — z 2 ) 2 , 1 ( 2z * 2 f 2 ). The 



asymptotic solutions for the periods are 

/1 = Xi(l + -xi +x 2 + -^x\ + ~^x x x 2 H ) 

/ 2 = Xxy/X^ 

4 23 

/3 = fi log(xi) + - + —x? + 3xix 2 H 

3 lb 

U = / 2 log(x 1 )+x 1 x|(- + ix 2 + ...) (A.77) 

5. C\ fl Coo. In the homogeneous coordinate (zi, z 2 , z 3 ) the divisor C\ is zi = 0, so 
the good local coordinates is (xi, x 2 ) = (z\, z 3 ). Since at the intersection z 2 7^ 0, 
we can choose z 2 = 1 and use the relation (jA.68|) to find (xi, x 2 ) = y). The 
asymptotic solutions for the periods are 

-3 23 1 
/1 = x x x 2 2 (1 - — Xi - -x 2 + • ■ ■) 

-1/ 1 67 2 7 1 2 

/2 = X 9 (1 X 2 Xi H XiX 2 x 9 + • • ■ I 

J2 2 v 4 144 1 24 32 2 y 

/3 = /llog(Xi)+X2 5 (x 2 -— X 2 + -XiX 2 --X2 + -..) 

/ 4 = (/1 - 2/2) log(x 2 ) + X 2 1 (~X 2 - ^X 2 + ^XiX 2 - -X^ + • • • ) 

b. In Cqo. In the homogeneous coordinates 5 1; z 2 and 5 3 , the divisor / is z 3 — 
2z,\ — 2z 2 = 0. At the intersection with / fl Coo the coordinate z 2 7^ 0, we can 
choose z 2 = 1 and use the relation ()A.b8|) to find the good local coordinates 

(xx, x 2 ) = (z 3 - 22, - 2z 2 , z 3 ) = ( i" 2 ^- 2 ^ 1) ( A . 78 ) 

z 2 z 2 
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The asymptotic solutions for the periods are 



ji = \Ai^ 2 i 1 + — ^ — ) 

. — I,,. 3 5 2 3 3 2 \ 

to = x 9 (1 H — Xi x, H iii 2 x H ) 

J 2 v 4 64 1 32 64 2 ; 

/ 3 = x 2 2 (xi + x 2 )(l + 3:1 - X<1 H ) 

J 4 = fl^Og{x 2 ) + y/x~ 1 X^(-^X l + ^X2 + ^xl + ^X 1 X2-^xl + ---) 

7. J fl Cqo. In the homogeneous coordinates zi, z 2 and 5 3 , the divisor J is — 
6Z1Z3 — 6,22-23 + 95 2 + \Az\z 2 + 9z% = 0- At the intersection with J fl the 
coordinate z 2 7^ 0, we can choose z 2 = 1 and use the relation ()A.68j) to find the 
good local coordinates 



(A. 79) 



(x 1 ,x 2 ) = (zl - Qziz 3 - Qz 2 z 3 + 9zl + 145i^ 2 + 9z 2 , z 3 ) 
1 - 6zi - 6z 2 + 9z\ + 14^1^2 + 9z| 1 

The asymptotic solutions for the periods are 

f - x\~H\ I 103 ~^ x I 8 + V *\ +...) 
Ji - x ± x 2 [L H — XiH — — x 2 H J 

/2 = ^ (1 + z^ II + z£z i ^ + ... ) 

64 16 

, 8192(43 + 13^),,, , v 
/s = (/1 2187 /2)log(x 2 ) 

-I, 256(7 + 4^) 1024(2 + 5^2i) 

+ X, X\ x 2 + • • • ) 

2 v 243 gl ; 

/4 = /llog , 1+ ,-i( 512 ' 23 - 10 ^> , 2 + ...) (A.™, 

A. 2 Analytic Continuation 

The periods of a Calabi-Yau manifold are integrals of the holomorphic three- form 
over the three-cycles. In the case of the Dijkgraaf-Vafa model, the integrals of the 
holomorphic three-form over the symplectic three-cycles reduce to integrals of a dif- 
ferential one-form over its branch cuts on the complex plane. For convenience we 
consider the cubic potential W(x) = \x 2 + |x 3 with the cubic coupling g set to one. 
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The A-cycle periods and B-cycle periods are 



Si 



So 



IT 



1 

2m 



X4 

dx\J [x — x\){x — x 2 ){x — x 3 )(x — X4) 



1 



2m 
1 
2~7ri 
1 

2~ni 



A 



3:3 
A 



fix a/ (x — — X2){x — X~sj(x — £4) 

dxyf {x — X\){x — X2){x — X$){x — X4) 

dxy/ (x — X\){x — x 2 ){x — x 3 )(x — x 4 ) 



(A.81) 



•»'2 



The asymptotic expansion of the periods around the origin (zi,^) = (0,0) was 
considered in |3|. Around this point the roots satisfy x\ < 22 < ^3 < £4 and the cuts 
between X\,x^ and between £3,0:4 shrink to zero sizes. It was found there that the 
asymptotic expansions of the periods are 



Si 

s 2 

2mlli 
2mU 2 



£1 
4 



-zi(2z 1 + 3z 2 ) + --- 

-^ + Iz 2 (3^ + 2^) + -.- 
4 o 

^4-M + 5- llog (£l)_fl 
3 2 &v 4 ; 4 

^-- + S 2 log(^' 
2 6 6V 4' 



A 3 



(Si + S 2 )\og(A 2 

Z-2 



+ 



S2) log(A^) 



(A.82) 



The Picard-Fuchs equations we derived can determine the constant term in the B 
cycle periods IT and the cut-off parameter is fixed to be Ao 
monodromy matrices around this point Zi — > z i: e 



|. We can find the 



2-1 



( 1 \ 
10 
10 10 

y 1 ) 



Mr. 



( 1 \ 
10 




\0 1 



1 
1/ 



(A.83) 



Since the integrals are done over symplectic cycles, the monodromy matrices are 
elements of the symplectic group Sp(4, Z) and satisfy the van Kampen relation 



M Z1 M Z2 



M Z2 M Z1 . 



Now we want to analytically continue the periods to other points in the complex 
structure moduli space. The analytic continuation will fix the symplectic basis of 
periods, which is not available by solving the Picard-Fuchs equation around these 
points. In order to do the analytic continuation, we must do the integrals in (JA.Slj) 
exactly. The A-cycle periods Si and the difference between the two B-cycle periods 



8 The periods are determined only up to a sign ambiguity due to the square root factors in the 
formula. Here we have taken the proper signs to match the convention of 0]. 
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Hi — n 2 can be written in terms of complete elliptic integrals of the first, second and 
third kinds, and one of the B-cycle periods involves incomplete elliptic integrals. 

We consider analytically continue the periods (jA.81|) to a singular point (zi, z 2 ) = 
(0, |) in the moduli space. This is the closest singular point to (-21,-22) = (0,0) in 
the moduli space. We will use the local coordinate (z%, z 2 ) around the intersection of 
the blow up divisor F\ : \ — z 2 = and divisor G\ : Z\ = as we did for solving the 
Picard-Fuchs equation 

We directly compute the asymptotic expansion of one B-cycle period n 2 and use the 
asymptotic expansion formulae of the complete elliptic integrals in ^1 to obtain the 
asymptotic formulae for other periods. 

For convenience we define a function in terms of complete elliptic integrals of the 
first kind K(k 2 ), the second kind E(k 2 ) and the third kind IT(a 2 , k 2 ) as the following 

2 i„2\ _ 1 



P(a\k 



48a 4 (l-a 2 ) 2 (a 2 -A; 2 ) 2 
x {a 2 [3a 8 + 3fc 4 - 4a 6 (l + k 2 ) - Aa 2 k 2 (l + k 2 ) 

+2a 4 (2 + k 2 + 2k 4 )}E(k 2 ) + [-3a 10 - 3A; 6 + a 8 (4 + hk 2 ) 
+a 4 k 2 (A - Wk 2 ) - 2a 6 (2 + k 4 ) + 3a 2 k A (l + 2k 2 )]K(k 2 ) 
+3[-a 12 - 5a 8 A; 2 + 5a A k 4 + k 6 + 2a 10 (l + k 2 ) 

-2a 2 k\l + k 2 )]Ii(a 2 ,k 2 )} (A.85) 

We start from the original matrix model point (zi,z 2 ) = (0,0) in the complex 
structure moduli space where X\ < x 2 < x 3 < x±. The expressions for the A-cycle 
periods can be found using formulae in jl^]. After some algebra we found 



Si — — dx\J (x A - x)(x - s 3 ) (a; - x 2 )(x - x\) 

27r Jx 3 



1 2(x 4 — x 3 ){x 2 — x 4 )(x 4 — ^l) 2 p ,%3 — (%4 — x 3 )(x 2 — X\) . 



2vr a/(x 4 - x 2 )(x 3 - Xi) x 3 - Xi (x 4 - x 2 )(x 3 - Xi) 

■y r x 2 

S 2 = — / dx\J (x 4 — x){x 3 — x){x 2 — x){x — Xi) 
2vr J Xl 

1 2(x 4 - x 2 )(x 2 - x 1 )(x 3 - x 2 ) 2 p X 2 - xi (x A — x 3 )(x 2 - x x ) 



2vr y/( X 4 - x 2 )(x 3 - Xi) x 3 - xx (x 4 - x 2 )(x 3 - Xi) 

(A.86) 
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and the difference between the two B-cycle periods is 

"X3 



27ri(n 2 — III) = / dx\/ (x 4 — x){x 3 — x)(x — x 2 ){x — Xi) 



X' 2 



_ 2(x 4 - x 2 )(x 3 - x 2 )(x 2 - x\) 2 p , x 3 - x 2 (x 3 - x 2 )(x 4 - x x ) 
a/(x 4 -x 2 )(x 3 - xi) x 3 - xi (x 4 - x 2 )(x 3 - xi) 

(A.87) 

We will take these exact formulae at the matrix model point and analytically continue 
to the local coordinate (jA.84|) . 

We can also directly compute the asymptotic expansion one of B-cycle periods Il 2 
around (z\, z 2 ) = (0, 0) as follows 



27iiU 2 = / dxJ{{x + I) 2 -- + z 2 ){x 2 -z 1 z%) 

Jy/Ph-I V 3 



^_ dxJ((x + iy-l + z 2 )(x- Z ^-^ + 0(z 4 )) 



(A.88) 



where I — y § — 25i^| + 25 2 . We can compute the integrals exactly for the first few 
leading terms written above, and expand around the cut off A = oo keeping only 
positive powers of A . 

We can now use the expressions for the periods (jA.86|) . (jA.87|) . (jA.88|) and obtain 
the asymptotic expansions to a first few orders 

* = £ A 

\/3 1 , 

S2 ~ -— fl + — h ~^M h 



1 >/3,, , , n , ^ 3) 
2vrffl 2 = ^ + ^- -^ + -^/ 3 -(S'i + 5 2 )log(12Ag) (A.89) 



n x = n 2 + — ^-(/ 4 + (3-log(2 4 3 3 ))/i) 

27T2 4 



A| A|_ J_ V3 
3 2 12 10£ 

where /i, f 2 , f 3 , / 4 are the asymptotic expansion of the solutions for Picard-Fuchs 
equation we found earlier 

1 _ 23 _ _ 5 
2"l6 Zl ~ 72^ ~ 46656 



fl — ZiZ 2 2 (l — 77— 7^1 — —ZiZ 2 — Z 2 + 



/ 2 = 5 2 2 (l + ^i+45 2 + ---) 

729 

h = 1 + 8H 2 3 - 1085i5 2 3 - — ~z\ + • • • 

o 

h = /ilog(5i) + 4(f -^2 + ---) (A.90) 
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These asymptotic expressions of periods are linear combinations of the 4 solutions 
to the Picard-Fuchs equations we found earlier, provided we choose the cut-off con- 
stant to be Ao = — |. Thus we have found the canonical basis for the symplectic 
cycles. It is easy to write down the monodromy matrices around this point 



M~ Z1 



( 1 








o\ 





1 








1 





1 





V° 








1 / 



/ -1 



M. 



Fi 



V 



2 







1 









-1 






2 

1 / 



(A.91) 



The monodromy around the singular divisor C\ is the same as before M Z1 = M Zl . We 
can also see that the monodromy matrices are elements of Sp(4, Z) group and satisfy 
the van Kampen relation .\//-. = Mp 1 M Zl . 

In general it is not easy to do the analytic continuation of periods. We use a 
numerical method to match the basis of solutions of Picard-Fuchs equation at different 
points of the moduli space. We consider the intersection of the singular divisor J and 
the blow up divisor i*\ : | — z 2 = 0. The local coordinate and the solutions for the 
Picard-Fuchs equation are 

Zi(6 - Uz 2 - 9z x ) 



Zl 



Uz 2 
9fe - i) 2 



Z-2 



- z 2 



(A.92) 



2 1/ 5 33 
h = z 1 z*(l-—z 1 + —z 2 + ---) 

_ 2/1 3. 43. , 
h = z 2 (l + -zx- — z 2 + ...) 

h = l- 648z 2 3 + --- 

h = f 1 log(z 1 ) + zl( 5 -^ + 32~z l -^~z 2 + 

15 7 



(A.93) 



Using numerical method we find the canonical basis of the periods as the following 

51 = (-0.14 -0.26i)/i + 0.082/4 

5 2 = (0.14 + 0.26i)/i + 2.6/ 2 - 0.048/ 3 -O.O82/4 
nx = -0.26z/i + 1.03i/ 2 - O.OI6/3 

IT 2 = 1.03z/ 2 - 0.016i/ 3 (A.94) 

The monodromy matrix of the divisor Fi is the same as we have derived in (jA.91|) . 
We can now write down the monodromy matrix of the singular divisor J by looking 
at the transformation around Z\ 



Mj 



-> z x e 2m 








I 1 


-2 


2 




1 


2 


-2 







1 







\ 





1 


/ 



(A.95) 



28 



For the singular divisor / : 1 — 2z\ — 2z 2 = 0, we find essential singularities instead 
of simple singularities. This can be seen from the asymptotic behavior of the solutions 
of the Picard-Fuchs equation at any point in the divisor. We find that the radius of 
convergence for the asymptotic expansion is zero, i.e. the series is always divergent. 
This is an interesting new feature of the Dijkgraaf-Vafa model. 



B Matrix model calculations 



In this Appendix we give some details of the matrix model calculations following the 
approach in |H1 Cj - The cubic matrix model can be expressed in the eigenvalues of the 
matrix 

$2 A 2 q\ 3 

W@) = tr (y + V } = + ~t ] (R96) 

t=i 

Then the partition functions Z and free energy F are 

Vol{U{N)) J N\(2tt) n J 1 1 1 K ' 

(B.97) 

where A(A) = rL<j(^* — ^j) ^ s ^ ne standard Verdermonde determinant from the 
measure of the matrix. We expand N\ eigenvalues around the critical points a\ = 
and N 2 = N — Ni eigenvalues around the critical points a 2 = — -. Suppose the 
fluctuation is m, z/j 

Xi = n h i = 1,2, •• -Ni 

X l+Nl = -- + Ui £ = 1,2,-- -JV 2 (B.98) 
9 

Then the potential and the Vandermonde determinant become 

N\ 3 N2 2 3 1 

wm = E(f + f- } - B| - 9 ~y )+ N > w (~b (R99) 

i=l * j=l -> " 



a 2 (a)= n ^i-^? n fa-"*) 2 n n ^-^+-t 

l<ii<i 2 <Ni l<h<h<N2 l<i<Ni l<j<N 2 " 



Now we can treat the expansion around this vacuum as a model with two ma- 
trices $1 with eigenvalues m and $2 with eigenvalues V{. The interaction terms 
rii<i<jVi rTi<j<jv 2 (A*^ — v i + g) 2 m ^ ne Vandermonde determinant can be exponen- 
tiated and written as potential for the two matrices, then the partition functions 
can be straightforwardly evaluated by expanding the potential and computing the 
expectations values of Gaussian matrix model [2]. We note the fluctuation around 

unstable critical point — has a wrong sign kinetic term — j-. However this model 



29 



is perturbatively well defined if we treat $x as a Hermitian matrix and analytically 
continue $2 to be a ant i- Hermitian matrix. Alternatively, one can also determine the 
perturbative part of the free energy by directly evaluating the Gaussian integral for 
various values of N\ and N 2 and solving for the coefficients in the perturbation series. 
Using this method we are able to push the computations of the free energy to th e 
eighth order, and provide many checks of the topological string calculations in (|4.(j7j) . 
The perturbative part of the free energy is the followings 



F pe rt = ~N 2 W{a 2 ) - 2N X N 2 \og{g) 
+ 



+ 



+ 



+ 



+ 



+ 



- 5Nx 2 N 2 + 5NxN 2 2 - In 2 3 ) + - ^ ) 



(^i 4 - y Ni " N2 + 59N ^ N ^ - y N,N 2 3 + ^N 2 4 ) + C-N, 2 ~ ^-N,N 2 + 7 -N 2 2 ) 
,56 «= 871 4 2636 , , 2636 2 , 871 , 56 * s 

( Y 1 ~ T 1 2 + "l - 1 2 ~ "1~ 1 2 + ~ NlN2 - y^ 2 ) 



9 4 



,332 Ar 3 923 Ar2Ar 923 ArAr2 332 Ar 3n .35 Ar 35 . 

+ ( Ni NiN 2 + NxN 2 N 2 ) + (—Nx N 2 ) 

9 3 3 9 J y 6 6 ' 

,512, = 6823 , T = 28765 , t4at 2 67310 , oo, 

N\ Ni 5 N 2 + Nx N 2 iVi 3 iV 2 3 ± ■ • • 

V 3 1 2 2 2 3 ' 



9 6 



+ ( iVi 4 Nx 3 N 2 + 153497V 1 2 7V 2 2 =f • • ■ ) + (3387V - 16327Vi7V 2 + 3387V r 2 



,1864 4 47083 

iVi — 

3 6 

(L^jVi 7 - 45118iVi 6 iV 2 + 2479807Vi 5 A^ 2 2 - 540378iVi 4 7V 2 3 ± ■ 
5 

54416 

+ ( Nx 5 - 1875287Vi 4 iV 2 + 5700667Vi 3 7V 2 2 t • • ■ ) 

5 

,66132 , , , ,5005 5005 , 

+ (— =— JVi 3 - 120880iVi 2 7V 2 ± ■ • • + - ~^N 2 ) 

5 3 3 



2 



9 W 



,65536. 1933906 A77Ar 13258178 Ar 6 A7 2 37761034 Ar 5 Ar 3 52780010 A7 4 Ar 4 
(— g— JV1 8 JVi 7 JV 2 + Nx 6 N 2 2 ^5^3 + JVi 4 jV a T • • • ) 

( 1762048^^ _ 129M560 VjVa + 5^776^ _ ^ W ± . . . ) 
( 1305280^4 _ 18059582^3^ + T . . . } 

( l^^ 2 _^ W2 + 1680m iV22) 

(59^92^ _ 680 9 87967 jVi8jV2 + 564130824^, _ ^^^^3 + , mmN ^ T . . . 

( Nx 7 - 97692942iVi 6 ^ 2 + 5373725407Vi 5 7V 2 2 - 1166263112iV 1 4 iV 2 3 ± • • • ) 

(12963696A^i 5 - 2444384277Vi 4 iY 2 + 745362156iVi 3 7V 2 2 qp . . . ) 
,86388296 . , 855302550 „, „ , , 8083075 . . 8083075 



g 12 



Ni N^N 2 ± • • • ) + ( Nx t, N 2 ) 

6 
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r, 11534336 Ar 10 1834216417 AT g Ar 5978643549 AT 8 Ar 2 „ i , r 7xr 3 

+ iV/ iVi 9 Ar 2 H iVi 8 7V 2 - 6444922816iVi 7 7V 2 3 

L v 3 12 4 " 

+ 147431576467Vi 6 A^ 2 4 - 19289163957iVi 5 7V 2 5 ± • • • ) 
(66841600^ - 4347551555^,^ + 29785674795^, 

-42151343305iV 1 5 ^ 2 3 + 58765399140iV 1 4 iV 2 4 T • • ■ ) 

(362264064*" - ^!^V* + 148144671957^,^, _ ^^3^3 ± . . . } 
^882324352^, _ 28221164683^ + 1853fl047 ^ lW T . . . j 
,693764720 2 3860943680 693764720 2x i lfi 

( g * 2 NiN 2 + iV 2 2 )J g 16 (B.101) 

This model also contains a non-perturbative part of free energy defined as the volume 
factor of the U(N) gauge group as in [27], where it was computed with the following 
result 

N 2 N 2 3 1 

Fn. P . = -flog(N 1 ) + ^\og(N 2 )--(N 2 1 +Nl)--log(N l N 2 ) 

00 B 1 1 

+2C / (-l) + V , 29 A ^^ + ^^) (B.102) 

This non-perturbative part of the matrix model has the correct universal leading 
behavior of Calabi-Yau near the conifold point of its moduli space, as first pointed 
out in in the context of c = 1 string compactified at self-dual radius. 
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